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Abstract. We report the discovery of multiple shells around the eruptive variable star V838 Mon. Two dust shells 
are seen in IRAS and MSX images, which themselves are situated in a shell of CO. This securely establishes 
V838 Mon as an evolved object. We revisit the light echo, which arises from scattering off the innermost resolved 
dust shell, to infer a distance to V838 Mon of <;5.5 kpc. The dynamical timescales of the ejected shells, location 
in the Milky Way and inferred luminosity are consistent with V838 Mon being a low-mass AGB star experiencing 
thermal pulses of which the 2002 eruption might have sent the star into the post-AGB phase. This scenario, 
however, is inconsistent with the presence of a (young and massive) B3 V companion. 
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1. Introduction 



V838 Monocerotis was discovered on January 6, 2002, 
(Brown 2002) as a new eruptive variable star, and peaked 
at my — 7 mag on February 5, 2002. It then declined in 
the optical whilst becoming brighter in the IR (Munari et 
al. 2002). This was accompanied by dramatic changes in 
the photosphere, making the object an L-type supergiant 
(Evans et al. 2003). The cause of the eruption of V838Mon 
and the nature of its progenitor are unclear. 

The distance to V838 Mon may be estimated from the 
evolution of the light echo, which was discovered in mid- 
February 2002 (Henden, Munari & Schwartz 2002). After 
first interpretations suggested a distance of d ~ 700 pc 
(Munari et al. 2002; Kimeswenger et al. 2002), subsequent 
revisions have led to ever increasing values for the distance 
of ~ 2.5 kpc (Wisniewski et al. 2003), kJ6 kpc (Bond et al. 
2003) and 8± 2 kpc (Tylenda 2004). Early estimates sup- 
ported the scenario proposed by Soker & Tylenda (2003) 
for the merger of two main-sequence stars, but later esti- 
mates allow for post-main-sequence scenarios. 

We present here evidence for multiple mass-loss events 
from the progenitor of V838Mon, which have occurred 
before the 2002 eruption. A dust shell is visible on IRAS 
images, and the central object is resolved by MSX. The 
double-shell system is situated within a larger shell which 
is visible on CO maps. We also revisit the evolution of the 
light echo and discuss the implications of our findings for 
the nature of the progenitor and eruption of V838 Mon. 



2. Fossil dust shells seen in infrared emission 

Warm dust, such as found in the vicinity of stars, emits at 
mid-IR wavelengths, whilst cold dust, such as encountered 
in the interstellar medium, emits at far-IR wavelengths. 
The sky around V838 Mon is included in both the IRAS 
all sky survey (12, 25, 60 and 100 /J.m), and the Galactic 
plane survey by the MSX in several mid-IR bands of which 
band A (8.3 /im) was the most sensitive. Both surveys were 
completed well before the recent outburst of V838Mon 
(1983 and 1996 for IRAS and MSX, respectively). 

2.1. IRAS data 

We retrieved IRAS images from the IRAS server 1 , of 
a 5° x 5° area centred on V838Mon (Equ J2000 RA= 
7 h 04 m 04 s .85, Dec= -3°50'51".l), corrected for zodiacal 
light. Because of the scanning, IRAS maps suffer from 
striping patterns. The worst of these was corrected by 
subtracting a 'flatfield', which was created by rotating the 
image through —7°, median-filtering within a box of 2' (12 
and 25 /jm) or 4' (60 and 100 /xm), averaging the rows, 
expanding the result to create a two-dimensional image 
and rotating through +7°. 

V838 Mon is clearly detected in all four IRAS bands as 
a compact source, surrounded by an elliptical ring of dif- 
fuse emission at ~ 7' from the central object in the North- 
South direction and ~ 10' East- West (Fig. 1). The image 
is an average of the four images, after scaling them to a 
similar dynamic range in intensity (this was done by sub- 
tracting the mean and dividing by the mode, which were 
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Fig. 1. IRAS composite image (see text) of a 50' x 50' 
field centred on V838Mon (cross). Apart from marginally- 
resolved emission within 1' of V838 Mon, the object is sur- 
rounded by a ring of diffuse emission (dotted ellipse). The 
positions of two bright stars are indicated for guidance. 



determined within an area encompassing the ring). There 
arc two bright IR objects nearby which are identified as 
the foreground star HD 296134 and the dusty carbon star 
CGCS 1538, but the diffuse IR source immediately North 
of HD 296134 does not have a known counterpart. 

The presence of a ring-like structure of emission sur- 
rounding the position of V838 Mon and its persistence in 
all four IRAS bands lend support to the significance of its 
detection. We developed a computer code to test this in 
a more objective manner. The procedure is as follows: (1) 
Along a circle with radius r, the median surface bright- 
ness and the difference between the 3 rd and 1 st quartiles of 
surface brightness values ("spread") are determined. The 
ratio of median and spread is a measure of the significance 
of the emission structure. Because a resolved structure is 
much more significant than the average signal-to-noise per 
pixel, the maps arc multiplied by a factor / = y/2nr/S, 
where 6 is the size of a resolution element for which we 
adopt a conservative 5 — 1'; (2) This is repeated taking, 
in turn, every pixel in the image as the centre of the cir- 
cle. Hence a map of the significance of ring-like emission 
around a certain position is produced; (3) This is repeated 
for increasing circle radii. 

This is a conservative approach, as the observed ring 
structure is clearly elliptical and a circular ring will there- 
fore never be able to capture all of the most significant 
structure. Yet it did recover a ring with a radius of 10' 
centred less than 1' North of the position of V838 Mon 
as the most likely (semi-)circular structure in the image, 
with a significance of 6.3 a. 
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Fig. 2. MSX band A image (see text) of a 25' x 25' field 
centred on V838 Mon, with an enlargement of the central 
3' x 3'. The position of the bright star HD 296134 is indi- 
cated for guidance. V838Mon is surrounded by emission 
within a little less than an arcminute from the central 
source. 

With an average surface brightness of S100 ~ 2 MJy 
sr _1 the ring around V838Mon is brightest at 100 /im, 
where the central source is quite weak. But with an av- 
erage surface brightness of S12 ~ 0.1 MJy sr _1 it is also 
clearly visible at 12 /im where the central source is much 
brighter than the ring; probably the dust in the ring is 
colder than the dust in the central source. The latter is 
marginally resolved at 12 and 25 fj,m. 

2.2. MSX data 

The MSX survey reaches a similar depth to the IRAS 
survey, but the spatial information within the MSX data 
is superior. The central source in the IRAS image is clearly 
resolved in the MSX data into several emission structures 
within a roughly circular area of ~ 1'.5 diameter (Fig. 2). 
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The band A (8.3 /im) image was corrected for striping 
along the Galactic longitude direction by subtracting an 
image created from an average of the columns in between 
the bright point sources seen in the picture. V838 Mon is 
located near the centroid of the emission complex. None 
of the emission knots near V838 Mon could be identified 
with stars in 2MASS JHK images (which do show what is 
believed to be the progenitor of V838 Mon) . 

The significance of the emission is estimated by rebin- 
ning the image such that the emission extends over only 
two (adjacent) pixels, each with a surface brightness of 
Sg.3 ~ 2.0 x 1CP 7 W m~ 2 sr _1 . Excluding the emission 
from the five obvious stars, the standard deviation of the 
fluctuations in the rebinned map is a ~ 5 x 1CP 8 W m~ 2 
sr _1 ; the emission associated with V838Mon is significant 
at the 1/2 x 4cr level, i.e. more than 6 a. Being the brightest 
square arcminute in the entire 25' x 25' field (apart from 
the five bright stars), and its disposition with respect to 
V838 Mon, lead us to conclude that it is very unlikely to be 
a chance superposition of fluctuations in the background. 
At nearly 4<r, the next-significant patch of emission in this 
part of the MSX band A image is the smudge at ~ 5' North 
of HD 296134, which may be identified with the bright IR 
source above HD 296134 in the IRAS image (Fig. 1). 

3. A large shell of carbon monoxide? 

Due to its proximity to the Galactic plane (I = 217°. 80, 
b = +1°.05) the sky around V838Mon was included in 
several surveys of emission from carbon monoxide (CO), 
of which the latest compilation is presented by Dame, 
Hartmann & Thaddeus (2001). Located in the direction 
just shy of what looks like a molecular cloud, V838 Mon 
seems to be situated off-centre within a bubble of CO emis- 
sion with a diameter of <~ 1° (Fig. 3). 

The eastern rim of the bubble ( "upper left" in Fig. 3) 
has only been covered in the Superbeam survey (Dame 
et al. 1987), at a relatively poor spatial resolution (A8 = 
0.5°). Its velocity-integrated flux level is ~ 0.1 K km s" 1 
per angular resolution clement. The western part adjoin- 
ing the molecular cloud was recorded at a higher spatial 
resolution (AO — 0.25°; May et al. 1993), and is also dis- 
cernible on 60 and 100 /im IRAS maps. The flux level of 
the edge of the molecular cloud delineating this part of the 
bubble is ~ 2 K km s _1 per angular resolution element. We 
note that in constructing the composite CO map, Dame 
et al. only integrated emission that exceeded 3x the rms 
noise level in the spectra, and hence any emission seen in 
their maps is statistically significant. 

The radial velocity of V838 Mon of v LSR ~ +50 km s" 1 
(Kipper et al. 2004) is rather high and suggests a location 
in the outskirts of the Milky Way disc (Fig. 3 in Dame et 
al. 2001). Given the fact that there is no other bubble to 
be seen within an area of at least 10 square degree, the 
chance for the line-of-sight towards V838 Mon to intersect 
the bubble (which covers an area less than 1 square degree) 
by sheer coincidence is certainly < 10%. New, sensitive 
high-spatial resolution CO (and radio continuum) maps 
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Fig. 3. A closeup of the Galactic CO map of Dame, 
Hartmann & Thaddeus (2001), centred on V838Mon 
which seems to be situated within a shell of CO emission 
<~ 1° in diameter. 

are desired in order to confirm (or otherwise) the physical 
association of the CO bubble with V838 Mon. 

4. The light echo revisited 

Since its discovery in February 2002, the light echo of 
V838 Mon's eruption has been monitored and used to in- 
fer the distance to V838Mon. The latest interpretation 
was presented by Tylenda (2004), who analysed images of 
the light echo taken in 2002. However, images are already 
available over a much larger timespan, and in light of our 
discovery of circumstellar dust shells around V838 Mon it 
is timely to revisit how the light echo can inform us about 
the distance and nature of V838Mon. 

We measured the diameter of the outer rim of the light 
echo visible on high quality images taken at the WHT, 
HST, SAAO and USNO over a period of more than two 
years since outburst (Table 1). There are many sources 
of uncertainty entering the error on the measured value, 
but the most important of these — irregularity and limited 
contrast of the rim — are small enough for a sensible anal- 
ysis, and the values given are accurate to ±1". Although 
the echo is clearly elliptical, the axial ratio (EW over NS) 
of 1.063 ± 0.006 is quite close to unity. 

Simple models for the evolution of the light echo, valid 
for the small angular extent of the echo, predict that the 
diameter grows with time since outburst, At, as: 



diameter = 2 arctan 



v / cAt(2r ± cAt) 
d 



(1) 



(see Fig. 4) where the minus sign is appropriate for a thin 
shell of radius r and the plus sign for a thin plane of 
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Table 1. Light echo measurements from images spanning 
the period March 2002 until February 2004. The times are 
with reference to t = 5 February 2002 (brightest visual 
maximum) . 



Observatory 


t — to diameter (") 
(days) East-West North-South 


ratio 


WHT 1 


53 


12.0 


11.2 


1.07 


HST 2 


86 


14.3 


14.0 


1.02 


SAAO 3 


101 


16.4 


15.3 


1.07 


HST 2 


106 


17.2 


16.2 


1.06 


HST 2 


211 


24.4 


23.1 


1.06 


SAAO 3 


224 


24.2 


23.1 


1.05 


HST 2 


267 


27.5 


25.3 


1.09 


USNO 4 


291 


29.3 


27.0 


1.09 


HST 2 


317 


29.7 


27.9 


1.06 


USNO 4 


625 


42.8 


40.2 


1.06 


HST 2 


735 


46.5 


43.6 


1.06 



^unari et al. 2002 



* http: / /hubblesite.org/ newscenter /newsdesk/archive / releases / 

2003/10 and 2004/10 

3 http://da.saao.ac.za/news/v838mon.jpg 

4 http://www. ghg.net/akelly/v8381ar2.jpg and v8381ar3.jpg 




Fig. 4. Geometry of the light echo effect when (1) scat- 
tered off a circumstellar shell with radius r (solid line) , or 
(2) scattered off a sheet of interstellar material at distance 
r in front of V838Mon (dashed). The time delay At with 
respect to the light traveling directly towards Earth (dot- 
ted) is related to the differential pathlcngth 5 = cAt. For 
the case d^>r (panel b) the pathlcngth (resp. 5\ or 82) has 
a simple dependence on r and projected distance R, with 
t&n6~R/d. 



scattering dust at a distance r in front of the light source 
and perpendicular to the line-of-sight (c is the speed of 
light and d is the distance to the light source). 

Both the optical images of the light echo and the in- 
frared emission as imaged by MSX show significant sub- 
structure in the scattering medium. Nevertheless, the evo- 
lution of the light echo diameter is remarkably monotonic 
in time (Fig. 5), which gives credence to the employment 



1 1 1 < 1 1 1 ^-^ 




200 400 600 

time [days] 



Fig. 5. Measured evolution of the light echo around 
V838 Mon, and the predictions for several simple models 
(Equation 1). 

of simple models to describe the global geometry of the 
scattering medium. 

At a distance d = 8 kpc, a small shell with radius 
r = 0.6 pc (projected radius 0/26) fits the early 2002 data 
but is inconsistent with the more recent data, whereas a 
larger shell with radius r = 0.9 pc (projected radius 0/39) 
fits the 2004 data but not the earlier data (Fig. 5). The 
only single shell solution is obtained by placing V838 Mon 
at d^S20 kpc, allowing for a larger shell (projected radius 
> 0/7). This would place V838Mon beyond the edge of 
the Milky Way disc, either within the Milky Way halo 
or within an as yet unrecognised dwarf galaxy. This is 
extremely unlikely, given the proximity of V838 Mon to 
the Galactic plane and its apparent association with a 
CO shell at the edge of a molecular cloud. 

Alternative solutions for d ~ 8 kpc include a hybrid 
solution of a thin shell (with radius r ~ 0.7 pc) and a thin 
plane, but the plane would have to be only r ~ 0.3 pc 
in front of V838 Mon and hence intersect the shell — an 
unlikely and rather ad hoc geometry. 

The data may be reconciled with a reasonable distance 
if we relax the assumption of a thin scattering medium. 
For a distance of d = 5.5 kpc, the most recent data may 
be reproduced by scattering under right angles (r = in 
Eq. 1) from dust at ~ 0.6' from V838Mon. Also at d = 5.5 
kpc, the early 2002 data may be reproduced by scattering 
off a shell with r ~ 0.3 pc (projected radius 0.19'). The 
rim of the light echo might first arise from scattering off 
the dense inner edge of the shell, but continue to advance 
through the shell outward after the reflection off the inner 
shell edge turns back into itself after about a year. The 
circumstellar material echoing the light of the February 
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2002 outburst is thus identified with the innermost re- 
solved dust shell as seen in the MSX image. For small 
dust grains (^50.1/im) scattering under angles greatly dif- 
ferent from 0° is inefficient; the outer rim of the echo will 
dim faster than by geometric dilution alone. 

This scenario will not work for a distance less than 
d <~ 5.5 kpc, because once scattering occurs under right 
angles we expect the projected expansion of the light echo 
to proceed linearly with time — and this is not yet seen. 
But it works for larger distances, where the scattering does 
not yet occur under right angles. We thus predict the light 
echo to enter a stage of linear expansion with time, which 
then yields an unambiguous value for the distance. 

5. The nature of V838 Monocerotis 

Despite a geometric method (the light echo) to estimate 
the distance to V838 Mon there remains uncertainty about 
its value as it depends on our interpretation of the circum- 
stellar dust environment and the scattering properties of 
the dust grains. Although the simplest model for the light 
echo would place V838Mon outside the Milky Way, we 
also present an interpretation where the echo proceeds 
through an extended shell, whereby the scattering angle 
approaches 90°. Accepting the latter as the more likely 
scenario, we adopt a distance of d = 5.5 kpc, but remind 
the reader that this is a lower limit. 

To determine the time at which the fossil shells were 
ejected we also need to know the expansion velocity. 
Although this might, in principle, be estimated for the 
inner shell from spectra of the scattered light echo, such 
an observation has yet to be made and will not be feasible 
for the outer dust shell. The CO shell expansion velocity 
might be measured from spectra of the CO line profile in 
the line-of-sight towards V838Mon, but this has not yet 
been feasible (cf. Rushton et al. 2003) due to the great 
amount of dilution rendering only the rim visible (where 
the radial velocity due to expansion is minimal) . 

We here calculate dynamical timescales by assuming a 
constant outflow velocity of v = 10 km s . This would 
yield 10 5 yr per projected parsec. A further complication 
in case of the CO shell is that V838 Mon is placed consid- 
erably off-centre and that the material might have slowed 
down as it had to plough through surrounding interstellar 
medium. For the subsequent dust shells this may not have 
been an issue as the CO shell would have swept clean a 
cavity for the dust shells to expand freely into, unless the 
removed interstellar medium was replenished by stellar 
mass loss in between eruptions. 

We thus estimate dynamical timescales, for the CO 
shell ti > 5 x 10 6 yr, for the outer dust ring t 2 ~ 1.4 x 10 6 
yr, and for the inner dust ring t 3 <~ 1.1 x 10 5 yr. If the 
current eruption is the fourth such event, then it seems 
that these events have occurred at an increasing rate. 

Placing V838Mon at a distance of (at least) 5.5 kpc 
makes its outburst peak at a bolometric luminosity L > 
10 5 L Q . The total energy radiated throughout the eruption 
(Bond et al. 2003) must then have been E > 10 38 J. The 



putative progenitor would have had L^IO 2 L Q . The 2002 
eruption and possibly earlier eruptions that produced the 
fossil shells might thus have been the result of helium- 
flashes (thermal pulses) on the AGB. Indeed, Kipper et al. 
(2004) report enhanced abundances of s-process elements, 
suggesting 3 rd -drcdgc-up on the thermal- pulsing AGB. 

Situated in the outskirts of the Galaxy, V838 Mon is 
likely to be metal poor (Kipper et al. (2004) estimate 
[Fe/H] ~ -0.4) and of low mass (M<>1 M Q ). Inter-pulse 
timescales are longer for lower-mass AGB stars, typically 
several 10 5 yr for a 1 M Q star, but decrease in time as the 
core mass grows (Vassiliadis & Wood 1993), in accordance 
with the above dynamical timescales of the shells around 
V838 Mon. An AGB star with little envelope mass will be 
severely affected by a thermal pulse at the base of its en- 
velope. We therefore propose that V838Mon experienced 
its fourth and final thermal pulse, resulting in the ejection 
of what may become a planetary nebula. 

Our scenario for a low-mass AGB star progenitor of 
V838Mon, however, is inconsistent with the presence of 
a B3V companion (Munari, Desidera & Henden 2002). 
Membership of a binary system in which the companion 
is a young and massive star implies that the erupting star 
must be young and massive too. Future studies must find 
a way of reconciling these two contrasting sets of observa- 
tional evidence. 
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